Efficient Leaping Methods for Stochastic Chemical Systems
نویسندگان
چکیده
Well stirred chemical reaction systems which involve small numbers of molecules for some species have a stochastic behavior and can be modeled by a continuous time, discrete state Markov process. An exact method for simulating the time evolution of the system is the Stochastic Simulation Algorithm, but this method is extremely slow for realistic biological systems. We propose an adaptive leaping method based on local error formulas and we compare the adaptive method with the exact method and with the explicit tau leaping method with constant time steps.
منابع مشابه
Slow Scale Tau-leaping Method.
For chemical systems involving both fast and slow scales, stiffness presents challenges for efficient stochastic simulation. Two different avenues have been explored to solve this problem. One is the slow-scale stochastic simulation (ssSSA) based on the stochastic partial equilibrium assumption. The other is the tau-leaping method. In this paper we propose a new algorithm, the slow-scale tau-le...
متن کامل$S$-Leaping: An adaptive, accelerated stochastic simulation algorithm, bridging $\tau$-leaping and $R$-leaping
We propose the S-leaping algorithm for the acceleration of Gillespie’s stochastic simulation algorithm that combines the advantages of the two main accelerated methods; the τ -leaping and R-leaping algorithms. These algorithms are known to be efficient under different conditions; the τ -leaping is efficient for non-stiff systems or systems with partial equilibrium, while the R-leaping performs ...
متن کاملThe numerical stability of leaping methods for stochastic simulation of chemically reacting systems.
Tau-leaping methods have recently been proposed for the acceleration of discrete stochastic simulation of chemically reacting systems. This paper considers the numerical stability of these methods. The concept of stochastic absolute stability is defined, discussed, and applied to the following leaping methods: the explicit tau, implicit tau, and trapezoidal tau.
متن کاملEfficient Finite-difference Methods for Sensitivity Analysis of Stiff Stochastic Discrete Models of Biochemical Systems
In the study of Systems Biology it is necessary to simulate cellular processes and chemical reactions that comprise biochemical systems. This is achieved through a range of mathematical modeling approaches. Standard methods use deterministic differential equations, but because many biological processes are inherently probabilistic, stochastic models must be used to capture the random fluctuatio...
متن کاملImplicit Simulation Methods for Stochastic Chemical Kinetics
In biochemical systems some of the chemical species are present with only small numbers of molecules. In this situation discrete and stochastic simulation approaches are more relevant than continuous and deterministic ones. The fundamental Gillespie’s stochastic simulation algorithm (SSA) accounts for every reaction event, which occurs with a probability determined by the configuration of the s...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2008